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Abstract. - Mechanical single molecule experiments probe the energy profile of biomolecules. 
We show that in the case of a profile with two minima (like folded/unfolded) periodic driving 
leads to a stochastic resonance-like phenomenon. We demonstrate that the analysis of such 
data can be used to extract four basic parameters of such a transition and discuss the statistical 
requirements of the data acquisition. As advantages of the proposed scheme, a polymeric linker 
is explicitly included and thermal fluctuations within each well need not to be resolved. 



Force spectroscopy has been used extensively to probe mechanically the interactions within 
a single biomolecule, like RNA and DNA, peptides or proteins by applying an externally 
controlled load on such a structure; for reviews see [1—7]. Since in these experiments typically 
the extension along one coordinate r is controlled, the data and results are rationalized in 
terms of the notion of a (free) energy profile G(r) along this coordinate. A major challenge 
for future improvements is to design these experiment is such a way that as much information 
as possible on such a profile can be extracted reliably from the data. 

In the paradigmatic situation of an experiment using an AFM cantilever, the biomolecule 
is attached via a polymeric linker to the tip of the AFM whose base is driven according to a 
certain protocol, see fig-d Up to now, in most applications in force spectroscopy, the force is 
applied according to a linear ramp protocol x(t) = xq + vt. Here, xq is the cantilever position 
at time t = and v is the constant ramp velocity. From the peak of the distribution of 
unfolding forces, one can extract two parameters of the underlying free energy landscape: the 
distance x u from the first minimum to the maximum and the spontaneous off-rate fcj 2 [8] , see 
fig- El These quantities characterize mainly the unfolding part of the profile. If one wants more 
information about the potential G(r), in particular on quantities determining the refolding 
events, one has to resort to a more sophisticated analysis. 

In the latest development, analyzing the distribution of the applied work spent in the 
unfolding process yields the full underlying profile using Jarzynski's relation [9-11]. In this 
context, we have shown that a periodic rather than a linear ramp improves the reconstruc- 
tion of the profile significantly [12]. Such an analysis based on Jarzynski's relation, however, 
requires first that the position of the tip z{t) is measured very accurately. Second, the re- 
constructed energy profile refers to the coordinate z which includes the linker molecule. It is 
nontrivial to extract from such a profile G(z) the quantity G(r) which is of main interest. 
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Fig. 1 Fig. 2 

Fig. 1 - Sketch of an atomic force spectroscopy experiment. The protein (thick lines) is connected via 
a polymeric linker (thin line) to the tip of an AFM cantilever at z(t) whose base at x(t) is externally 
controlled. 

Fig. 2 - Schematic view of the underlying free energy potential G(r) as a function of the reaction 
coordinate r. The first minimum represents the folded state, whereas the second shallow minimum 
represents the unfolded state of the biopolymer. :r u and xt denote the respective distances from the 
minima to the barrier and k°j the spontaneous transition rates from state i to j. 



In this letter, we propose a method on an intermediate level of sophistication, which yields, 
apart from x u and fe° 2 , the corresponding two quantities Xf and fc 21 of the refolding transition. 
We employ a periodic ramp which leads to a stochastic resonance-like phenomenon [13]. As 
advantages of our new scheme, we include the linker and do not require to analyze the full 
stochastic trajectory z(t) as it is necessary in approaches based on Jarzynski's relation. It is 
sufficient to observe the sequence of folding and unfolding events reflected in this coordinate 
as explained below. Our proposal extends a previous approach where spontaneous transitions 
under constant force are used to extract these parameters [14]. The main advantage of our 
dynamical approach is that stochastic resonance leads to an enhancement of the transition. 
Moreover more control parameters such as amplitude and frequency of the driving are available 
to optimize the data acquisition. Finally no feed-back loop in order to impose a constant force 
is necessary. 

We first recall that the spontaneous transition rates k® 2 and fc 21 for unfolding and refolding 
are modified by a force F according to [8, 15, 16] 

fci2 = ^12 exp[/3x u F(z, L)\ (1) 
fcai - k° 21 cxp[-f3x f F(z,L)} (2) 

where (3 is the inverse temperature. The force F(z,L) is transmitted from the cantilever to 
the biomolecule by a polymeric linker which we model by a worm-like chain as [17] 

f ^ l) = wX a { i- z ,l ? -\ +Z l) < 

where L p is the persistence length. The contour length L depends on the configuration. In 
the folded state, it is L — L , in the unfolded state L = L + AL, where AL is the length 
gain of the contour length upon unfolding. 



(3) 
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In each state, the force F c measured by the cantilever with harmonic spring constant k c 
has to balance the polymeric force 

F c = k c (x-z) = F(z,L) . (4) 

Thus, the z-coordinate can be calculated for given externally controlled cantilever base position 
x. In general eq. (@J can not be solved analytically for z{x). In force spectroscopy on 
biopolymers, however, the linker is almost stretched out for a typical unfolding force of about 
50-150 pN. Eqs. © and (@J then reduce to the cubic equation 

k c (x-L + 6) = , (5) 

for S = L — z, which denotes the length stored in the fluctuations. Hence the functions 
z = z{x, L) or (5 = 5(x, L) also depend on the current length of the polymer. 
We now apply a periodic protocol of the cantilever position x according to 

x(t) = xq + i a cos(cjt) , (6) 

with the offset extension xq, the amplitude x a and frequency to [12]. In general, the periodicity 
of x(t) leads to an anharmonic periodicity of S(t) via eq. JSJ, but the leading term for x & <C a?o 
is given by the first harmonics as 

8{t) w S +<5 a cos(wt) , (7) 
where Sq is the solution of eq. JSJ with x = xq. The amplitude <5 a is given by 

Sa = ~3 + 2(x a -L)/S ■ (8) 

To first order in 5 a /5o, the force acting on the biomolecule exhibits the same periodicity 

L 2 ( 28 \ 
F(8, L) » ^1 - cos M )j = F + F a cos(orf). (9) 

For the functions Fq and F a we will have to distinguish between the folded the unfolded 
configuration. 

The probability $(t) that the molecule is in the folded state obeys the master equation 

^$> = -(k 12 + k 21 )$ + k 2 i, (10) 

where the rates depend on the force through and ■ For small loading, an analytical 
solution can be calculated in linear order in x^/xq by splitting the solution of the master 
equation (|Trj|) into a time-independent part $o and a time-dependent part II(t) as 

$(t) = $ + n(*) . (ii) 



The stationary part is given by 
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with 

«12 = k° 2 exp[/3a; u Fo((5o, L )] and k 2 i = k 21 exp[-/3xfF (5o, L + AL)] , (13) 

which are the time-independent contributions to the transition rates due to the preloading. 
On the intrinsic time scale 1/A = (k 12 + k 2 i) _1 , the time dependent part of $(i) approaches 
a limit cycle 

n(t) ee n cos(wt — A) = a cos(wt - A) , (14) 

with 

a ee -(/3z u F a (L ) + 0XiF a {L o + AL))k 12 $ (15) 

and the phase shift 

A ee arctan(w/ A) . (16) 

The solution of eq. (|1(J|> is thus a periodic response of the probability of occupation <I> of the 
folded state. 

For fixed xo and fixed u> the solution lllfl yields three parameters a, A and $o for the 
four unknowns x u , Xf , k® 2 , k 21 . F° r a determination of these crucial parameters, more data are 
needed, which can easily be obtained by varying the offset extension Xq and the frequency to. 
Such a variation leads to an interesting resonance phenomenon. 

It is one of the well known features of periodically driven stochastic transitions, that 
in contrast to a deterministic mechanical system, the response amplitude depends on the 
frequency in a non-resonant manner [13]. For low frequencies, where the stochastic transition 
can follow the externally driven loading, the amplitude is larger than for large frequencies, see 
eq. I|14|) . Both, the response amplitude Ilo of the periodic part of the occupation probability 
and the phase shift A show a clear maximum as a function of the offset xo, see figs. [3] and 0] 
The maximal phase shift A* follows from the condition k' 12 + 1' 12 = 0, where the prime denotes 
the derivative with respect to xq . This yields the ratio of the geometric parameters X[ and x n 
as 

Xf_ = Ki2^(Xo) 

x u K 2 iF^(L +AL) ' 1 ' 

where Fq(Lq) and Fq(Lq + AL) are the derivatives of the force with respect to xq with folded 
contour length and unfolded contour length respectively. With eq. 101, we get explicitly 

L 2 

FL(L) = 5 , (18) 

2PL p 5* 2 (6* + 2(S*+x*-L)) 

where Xq is the position of the maximum of the phase shift and Sq — So(xq, L). Eq. (|17|l then 
yields the ratio x u /xf of the underlying free energy profile as a fourth relation which can be 
used to determine x u ,Xf, kf 2l k 21 . Of course, many other fitting schemes can be envisaged to 
extract these crucial parameters. 

Up to now, we have not yet addressed explicitly the question of how to obtain the proba- 
bility $(£) from real measurements and how these experiments have to be designed to deliver 
the necessary data most efficiently. In principle, the occupation probability $(£) has to be 
measured via an ensemble average of trajectories. The ensemble average, however, can be 
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Fig. 3 - The analytical solution of the response amplitude of the occupation probability of the folded 
state Ho as a function the offset extension xq. The parameters are k^i = 3s _1 , fc? 2 = 4.6 • 10 _4 s _1 , 
x u — 0.3 nm, Xf = 1 nm, L p — 0.42 nm , Lo = 40 nm , AL = 28.5 nm, u = Is - , x a — 0.5 nm, 
fc c = 5pN/nm, taken from [18] for a tenascin FN-III segment. 

Fig. 4 - The phase shift A of the occupation probability over the offset extension xo for the same 
parameters as in fig-B] 



replaced by a phase sensitive time average of a single trajectory, because <&(t) is periodic with 
period time T — 2it/lu after a relaxation time for the offset expansion of the order of t ~ 1/A. 
Comparing these two averaging schemes, the advantage of the phase sensitive time average is 
obviously that the relaxation time has to be accounted for only once. 

According to our assumptions, the trajectories z(t) are bistable, so the conformation 
(folded or unfolded) of the biopolymer can be determined unambiguously from z(t) since 
each transition leads to a jump in z(t) which is larger than the thermal fluctuations within 
each minimum. We assign a weight p(t) = 1 if the molecule at time t is in the folded state 
and p(t) = if it is in the unfolded state. The probability of occupation $(t) of the folded 
state can then be estimated by 

= , (19) 

where N = T tot /T is the number of periods with the total measuring time T tot and p n (t) = 
p(t + nT) and < t < T. We now optimize uj and T tot for the least statistical error. 
The autocorrelation of p(t) up to first order in x&/xq is given by [13] 

cov(p(t),p(t + r)) ^ («p(t)) - p{t)){{p{t + r)> - p{t + r))) « , (20) 

with the single measurement error of a two state system a^(t) = <1> (t ) ( 1 — $(£)). We consider 
the A^-dependent stochastic variable $^(i), which fluctuates around the mean average value 
$(t) with a variance o% N = J2n=o ^pi^N / dp n ) 2 +Y<n^m( d ® n /dp n )(d$ N /dp m )cov{p n ,p m ), 
where we omitted the i-dependence. Together with eq. I|20() we get approximately a geometric 
sum, which yields 
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Fig. 5 - The minimum total measuring time T m i n over the offset extension xq for the optimized 
frequency u>* ~ A and the same parameters as in figs. |3 and 0] 

The square root of the variance of this stochastic variable can now be identified with the 
measuring error, which should be smaller than the response amplitude ITo. Thus, we get the 

condition \ cr% < |IIo|. The total measuring time must therefore obey 



where we have replaced the phase sensitive $(t)(l — by <E>o(l — $o)> which is permissible 
in zeroth order of O(x a /xo). The right hand side of equation (|22[l has a minimum at lo* ~ A, if 
the exponential correction in the denominator is neglected. The total measuring time for the 
optimized frequency lu* is T m i n = Ttot(u>*) = 47tKi2K2i/(Aa 2 ), which strongly depends on the 
offset extension xo, see fig. El The minimum around Xq is a signature of a stochastic resonance- 
like phenomenon. Usually in stochastic resonance, the potential and a periodic driving input 
is given while the noise strength is controlled. One then finds an optimal noise strength for 
the best amplification of the input signal. In our case the strength of the noise is fixed by 
the temperature, whereas the driving force can easily be manipulated by adjusting offset Xq, 
amplitude x a and frequency to. Figs. and [S] show that where the response amplitude n as 
a function of offset x$ is the largest the total measuring time T min is the smallest. Hence, the 
signal to noise ratio is the best for this choice of offset and optimized frequency. Moreover, 
it should be noted that typically in stochastic resonance the bistable potential is symmetric 
whereas in our case it is intrinsically asymmetric due to the molecular interactions. 

In summary, we have shown that the two spontaneous transition rates k® 2 , ^21 an d the 
two most relevant length scales x u and Xf of a biomolecular conformational transition can be 
extracted from data obtained by periodically driving these transitions. Our analysis focuses 
on the measured occupation probability of one state, which exhibits a stochastic resonance-like 
phenomenon as a function of the control parameters of a typical AFM experiment. Compared 
to more sophisticated approaches based on Jarzynski's relation, our scheme does not require to 
resolve thermal fluctuation within each minimum. Moreover, the linker molecule is explicitly 
included. Future applications of our scheme will include more complex transitions and linear 
sequences of similar transitions in the un- and refolding of multi-domain biopolymers. Most 
important, however, are experimental tests of periodic loading which we hope to encourage 
with this study. 




2tt$ (1 - $ )(^ 2 + A 2 
a 2 uj{l - e -A27r/^) 



(22) 
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